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In general, the power-law degree distribution has profound influence on various dynamical pro- 
cesses defined on scale-free networks. In this paper, we will show that power-law degree distribution 
alone does not suffice to characterize the behavior of trapping problem on scale-free networks, which 
is an integral major theme of interest for random walks in the presence of an immobile perfect ab- 
sorber. In order to achieve this goal, we study random walks on a family of one-parameter (denoted 
by q) scale-free networks with identical degree sequence for the full range of parameter q, in which 
a trap is located at a fixed site. We obtain analytically or numerically the mean first-passage time 
(MFPT) for the trapping issue. In the limit of large network order (number of nodes), for the whole 
class of networks, the MFPT increases asymptotically as a power-law function of network order with 
the exponent obviously different for different parameter q, which suggests that power-law degree dis- 
tribution itself is not sufficient to characterize the scaling behavior of MFPT for random walks, at 
least trapping problem, performed on scale-free networks. 

PACS numbers: 05.40.Fb, 89. 75. He, 05.60.Cd, 89. 75. Da 



I. INTRODUCTION 



As a fundamental stochastic process, random walks 



of trapping is determining its efficiency, namely, showing 
the dependence relation of MFPT on the size of the sys- 
tem where the random walks are performed. Previous 
studies have provided the answers to the corresponding 
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In this paper, we study the trapping problem on the 
a class of scale-free networks with the same degree se- 
quence, which are dominated by a tunable parameter 
q [7 ]. We determine separately the explicit formulas 
of the mean first-passage time for the two limiting cases 
of q = 1 and q = 0. Wc show that in both cases the 
MFPT increases as a power-law function of the network 
order, with the exponent less than 1 for q = 1 and equal 
to 1 for q = 0. We also study numerically the MFPT for 
the case of < 5 < 1, finding that it is also a power- 
law function of network order with the exponent 6{q) 
depending on parameter q. We demonstrate that in the 
full range of < g < 1, 0(g) is a decreasing function of 
q, which belongs to the interval [|^, l] . Our findings in- 
dicate that the power-law degree distribution by itself is 
not sufficient to characterize the trapping process taking 
place on scale-free networks. 
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FIG. 2; (Color online) Illustration of the first several itera- 
tions of the network for a particular case q — 1. 



II. THE SCALE-FREE NETWORKS WITH 
IDENTICAL DEGREE SEQUENCE 




FIG. 1: Iterative method of the network construction. Each 
edge is replaced by either of the connected clusters on the 
right-hand side of arrows with a certain probability, where 
black squares represent new nodes. 

The networks in question are built iteratively [? ], 
see Fig. ??. We represent by Hn (n > 0) the networks 
after n iterations (the number of iterations is also called 
generation hereafter) . Then the networks are constructed 
as follows. For n — 0, the initial network Hq consists of 
two nodes connected to each other by an edge (a link). 
For n > 1, Hn is obtained from Hn-i- That is to say, 
to obtain one can replace each link existing in Hn-i 
either by a connected cluster of links on the top right of 
Fig. ?? with probability q, or by the connected cluster 
on the bottom right with complementary probability 1 — 
q. Repeat the growth process n times, with the graphs 
obtained in the limit n — > 00. In Figs. ?? and ??, we 
present the growing processes of two special networks 
corresponding to q — and q — 1, respectively. 

Let L(ri) be the number of nodes created at generation 
n, and the total number of all edges present at gen- 
eration n. By construction, we have En = 4 Con- 
sidering the initial condition _Eo = 2, it leads to En — 4". 



Since each existing edge at a given generation will create 
two new nodes in at the next generation, then, at each 
generation rii {rii > 1) the number of newly introduced 
nodes is L{ni) — 2£'„._i — 2x4"'~^. Thus, at generation 
n the network order is 

" 9 
K= 5]L(n,) = -(4" + 2). (1) 

ni=0 

Let ki{n) be the degree of a node i at generation n, 
which was created at generation ni {rii > 0). Then, 

fc,(n) = 2"-"'+i. (2) 

From Eq. (??), it is obvious that after each new iteration 
the degree of a node doubles, i.e., 

h{n)^2h{n-l). (3) 




n=2 n=3 



FIG. 3: Sketch of the iteration process of the network for the 
limiting case of q = 0. 

The networks considered exhibit some interesting topo- 
logical properties. Their nodes have same degree se- 
quence (thus the same degree distribution), independent 
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of the value of parameter q. Concretely, the networks 
have a power-law degree distribution with the exponent 
7 = 3 [? ]. On the other hand, since there is no tri- 
angle in the whole class of the networks, the clustering 
coefficient is zero. Although the degree distribution and 
clustering coefficient do not depend on the parameter q, 
other structural characteristics are closely related to q. 
For example, ior q — 1, the network is reduced to the 
(2, 2)-flower introduced in [? ]. In this case, it is a 
small world, its average path length (APL), defined as 
the mean of shortest distances between all pairs of nodes, 
grows logarithmically with the network order [? ] ; at the 
same time, it is a non- fractal network [? ? ]. While for 
g = 0, it is exactly the hierarchical lattice that was pro- 
posed by Berker and Ostlund [? ] and was extensively 
studied by many authors [??????]. For this case, 
the network is not small-world with the APL increasing 
as a square power of the network order [? ? ] ; moreover, 
it is fractal with the fractal dimension rfs = 2. When 
q increases from to 1, the networks undergo a tran- 
sition from fractal to non-fractal scalings, and exhibit a 
crossover from 'large' to small worlds at the same time [? 
]; these similar phenomena are also observed in a family 
of treelike networks [? ]. 

The peculiar topological features make the networks 
unique within the category of scale-free networks, since 
these particular structures strongly affect the dynamical 
processes defined on the networks. For instance, different 
thresholds of bond percolation were recently observed in 
the networks, which implies that power-law degree dis- 
tribution alone does not suffice to characterize the perco- 
lation threshold on scale-free networks under bond per- 
colation [? ? ]. In what follows, we will study random 
walks with a single immobile trap on the networks. We 
will show that the degree distribution is not sufficient 
to determine the scalings for MFPT of trapping process 
occurring on the networks under consideration. 



III. RANDOM WALKS WITH A FIXED TRAP 

In this section, we study the so-called simple discrete- 
time random walks of a particle on network At each 
time step, the particle (walker) jumps from its current 
location to one of its neighbors with equal probability. 
In particular, we focus on the trapping problem, i.e., a 
special issue for random walks with a trap positioned at 
a given node. To this end, we first we distinguish differ- 
ent nodes in by labeling them in the following way. 
The two nodes in Hq have labels 1 and 2. For each new 
generation, we only label the new nodes created at this 
generation, while we keep the labels of all pre-existing 
nodes unchanged. In other words, we label sequentially 
new nodes as M + 1, Af + 2, . . . , M + AM, where M is 
the number of the old nodes and AM the number of 
newly-created nodes. In this way, every node is labeled 
by a unique integer, at generation n all nodes are labeled 
from 1 to = 1(4" -|- 2). Figures ?? and ?? show how 
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FIG. 4: labels of all nodes of H3 in the case of g = 1. 
1 




FIG. 5: labels of all nodes of H3 for the particular case of 
q = 0. 



the nodes are labeled for two special cases of g = 1 and 
q = 0. 

We place the trap at node 1, denoted by ix- At each 
time step, the particle, starting from any node except the 
trap ir, moves uniformly to any of its nearest neighbors. 
It should be mentioned that, due to the symmetry, the 
trap can be also situated at nodes 2, 3, or 4, which has not 
any effect on MFPT. The special selection we made for 
the trap allows to address the issue conveniently. Partic- 
ularly, this makes it possible to analytically compute the 
MFPT for the two deterministic networks corresponding 
to q = 1 and q = (details will be discussed below), be- 
cause of their special structures and the convenience of 
identifying the trap iy since the first generation. 

As mentioned above, one of the most important quan- 
tity characterizing such a trapping problem is the FPT 
defined as the expected time a walker takes, starting from 
a source node, to first reach the trap node. The signif- 
icance firstly originates from the fact that the first en- 
counter properties are relevant to those in a plethora of 
real situations [? ], including transport, disease spread- 
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ing, target search, and so on. On the other hand, many 
other quantities can be expressed in terms of FPTs, and 
more information about the dynamics of random walks 
can be extracted from the analysis of FPTs [? ] . Finally, 
the average of first-passage times, i.e., mean first-passage 
time (MFPT), measures the efficiency of the trapping 
process: the smaller the MFPT, the higher the efficiency, 
and vice versa. In the following, we will determine the ex- 
act solutions to MFPT for some limiting cases, as well as 
the dependence relation of MFPT on the network order. 

Let t/"' be the FPT for a walker initially placed at 
node i to first reach the trap ix in Hn. This quantity can 
be expressed in terms of mean residence time (MRT) [? 
? ] , which is defined to be the mean time that a random 
walker spends at a given node prior to being absorbed 
by the trap. Actually, the MRT is the mean number of 
visitations of a given node by the walker before trapping 
occurs. 

It is known that the trapping problem studied can be 
described by a Markov chain [? ], whose fundamental 
matrix is the inverse of matrix B„ that is a variant of 
the normalized Laplacian matrix [? ] L„ for H„, and can 
be obtained from L„ with all entries in the first row and 
column (corresponding to the trap node) setting to ze- 
ros. The entry {b~^)ij of the fundamental matrix (B„)~^ 
expresses the mean number of visitations of node j by 
the particle, starting from node i, before it is eventually 
trapped. Thus, we have 

= f:{bn% ■ (4) 

Then, the mean first-passage time, (T')„, which is the av- 
erage of t/"'' over all initial nodes distributed uniformly 
over nodes in iJ„ other than the trap, is given by 

(^)« = v^E^^/''-^EE(^n%- (5) 

1=2 1=2 J=2 

Equation (??) shows that the problem of determining 
(r)„ is reduced to computing the sum of all elements 
of the fundamental matrix (B„)~^. Although the ex- 
pression of Eq. (??) seems compact, the complexity of 
inverting L„ is 0{V,^). Since the network order increases 
exponentially with n, Eq. (??) becomes intractable for 
large n. Thus, restricted by time and computer mem- 
ory, one can obtain (T)„ through direct calculation from 
Eq. (??) only for the first iterations. It would be satisfac- 
tory if good alternative computation methods could be 
proposed to get around this problem. What is encour- 
aging is that the particular construction of the networks 
and the special choice of the trap location allow to calcu- 
late analytically MFPT to obtain a closed-form formula, 
at least for the two special cases of g = 1 and q = 0. 
The computation details will be provided in the follow- 
ing text. 



A. Case of g = 1 

We first establish the scaling relation governing the 
evolution for T^^"^ with generation n. In Table ??, we 

list the numerical values of T!^"-* for some nodes up to 
n = 6. From the numerical values, we can observe that 
for a given node i, the relation T"/""*"^' = 3T/"' holds. 
That is to say, upon growth of the network from gener- 
ation n to generation n -I- 1, the trapping time to first 
arrive at the trap increases by a factor 3. This is a basic 
characteristic of random walks on when g = 1 , which 
can be established from the arguments below [???]. 

Consider an arbitrary node i in _ff„ of the q = I case, 
after n generation evolution. From Eq. (??), we know 
that upon growth of the network to generation n + 1, the 
degree, ki, of node i doubles, namely, it increases from 
ki to 2ki. Among these 2 fcj neighbors, one half are old 
neighbors, while the other half are new nodes created 
at generation n+1, each of which has two connections, 
attached to node i and another simultaneously emerging 
new node. We now examine the standard random walk 
in Let X be the FPT for a particle going from 

node i to any of its ki old neighbors; let Y be the FPT 
for going from any of the kj new neighbors of i to one 
of the ki old neighbors; and let Z represent the FPT for 
starting from any of new neighbors (added to the network 
at generation n+1) of an old neighbor of i to this old 
neighbor. Then we can establish the following backward 
equations: 

(X= l + i(i + y), 

I Y = i(l + X) + i(l + Z), (6) 

[z= i + i(i + y). 

Equation (??) has a solution X ^ 3. Thus, upon the 
growth of the network from generation n to generation 
n+1, the first-passage time from any node i to any node 
j (both i and j belong to _ff„) increases by a factor of 3. 

That is to say, = 3i;^"\ which will be useful for 

the derivation of the exact formula for the MFPT below. 

After obtaining the scaling of first-passage time for old 
nodes, we now derive the analytical rigorous expression 
for the MFPT (T)„. Before proceeding further, we first 
introduce the notations that will be used in the rest of 
this section. Let A„ denote the set of nodes in and 
let An stand for the set of those nodes entering the net- 
work at generation n. For the convenience of computa- 
tion, we define the following quantities for 1 < m < n: 

t(") _ \ " rp(n) ,„\ 

-'TO,t0t — 2-^ -'j ' \' ) 

and 

By definition, it follows that A„ = A„ U A„_i. Thus, 
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